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Abstract 



Leptonic decay widths and leptonic decay constants of light vector mesons 
and weak leptonic decay widths and weak decay constants of light and heavy 
■ pseudoscalar mesons have been studied in a field- theoretic framework based on 

£SJ , the independent quark model with a scalar- vector power-law potential. The 

results are in very good agreement with the experimental data. 
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1 Introduction 

Though quantum chromodynamics is considered to be the underlying theory of 
D . strong interaction between quarks and gluons at the structural level of hadrons, 

many low-energy phenomena such as spectroscopy, static electromagnetic properties, 

•rH , 

can not be explained by first-principles application of QCD. Therefore one needs to 

5-H ' 



resort to phenomenological models. For example, a potential model with an equally 
mixed scalar- vector harmonic potential 1 of independent quarks in a relativistic Dirac 
framework has been used to study several low-energy phenomena in the baryonic 
sector such as octet baryon masses 2 , magnetic moments 3 , weak electric form factors 4 , 
nucleon electromagnetic form factors and charge radii 5 . In addition to the harmonic 
potential, with power-law potential such as r v {y > 0) heavy- heavy and light-heavy 
quarkonium states have been studied and the results are in good agreement with 
experimental results 6 . This model has also been successful in explaining pion mass, 
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its decay constant 7 and the radiative decay 8 of ordinary light and heavy mesons. 
Because of this wide range of applicability of the model to both baryons and mesons, 
it has proved to be a rather simple and successful alternative to the cloudy bag 
model 9 . The purpose of this work is to extend its applicability to the study of 
the leptonic decay of vector mesons such as p, u, 4> and weak leptonic decay of 
light and heavy pseudoscalar mesons for an equally mixed scalar-vector power-law 
potential such as Ar u . The leptonic decay width of heavier vector mesons in the 
charm and bottom quark sector have been studied in the nonrelativistic approach 
through the Van Royen-Weisskopf formula with radiative corrections 10 . But the 
same approach is not suitable for ordinary vector mesons in the light flavor sector, 
where the constituent quark dynamics is more relativistic. On the other hand, 
for weak leptonic decays, while many nonrelativistic quark model calculations 11 
suggest that fx > Id > Ib {Im is the weak leptonic decay constant), some of 
the models based on QCD sum rules 12 and lattice calculations 13 predict more or 
less a constant /m between K and B mesons. Capstick and Godfrey 14 calculate 
the hadronic matrix elements for the relativitized quark model expression for f\j- 
They find f Bc > Jd s > Ik > I'd > Jb s > I'b > U, but the calculated value of 
the ratio fx/ fir — 1-75 is much higher than the experimental value of 1.22. The 
leptonic decay widths and decay constants fy of the light vector mesons have been 
calculated by using an equally mixed scalar vector harmonic potential 15 . Also, weak 
leptonic decay constants, Jm, of pseudoscalar mesons have been calculated using the 
same potential and are found to satisfy /# c > fr> s > /d > fx > /b s > fw > /s 16 - 
Potential used in Ref. 15 and 16 is harmonic. Since the potential between quarks is 
weaker than the harmonic potential, the potential used in this study is Ar - 2 + Vq 
which gives very good results for heavy-heavy and light-heavy quarkonium states 
and radiative decay of ordinary light and heavy mesons 6 ' 8 . Both in Ref. 15,16 and 
in this study, the potential does not include Coulomb term because Coulomb-like 
vector interaction are believed to have less prominent role for the light mesons. 
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2 Potential model 



The quark-confining interaction in a hadron, which is believed to be generated by 
the nonperturbative multigluon mechanism, is not possible to calculate theoreti- 
cally from first-principles within QCD. Therefore, from a phenomenological point of 
view, the present model assumes that the quark and antiquark in a hadron core are 
independently confined by an average flavor-independent potential V(r) of the form 



V(r) = ^(l + P)(Ar u + Vq), A>0 and v > 0. 



For this potential, Dirac equation can be written as 



[a-p + [3m + V(r)] *(r) = EV(r) 



(2) 



where a and (3 are Dirac matrixes. Eq. (2) has two solutions with positive and 
negative energy given respectively in the forms 



ig\(r)/r 
o ■ rf A (r)/r 

o ■ rf A (r)/r 
g\{r)/r 



17a (r) 



U(r) 



(3) 



(4) 



where A = (nljm) represents the set of Dirac quantum numbers specifying the 
eigenmodes. The spin angular parts U\(r) and U\(r) are described as 



U ljm (r) = Yl C'm^l^r^. 



TO; ,TO S 



U ljm (r) = (-iy +m ' l U^ m (r 



(5) 
(6) 



Substituting Eq. (3) or Eq. (4) into Eq. (2), one obtains (for n = 0, I = 0) 

g A {r) = X q (E-m-V )gA(r), (7) 



2 dr 2 



— — + X q Ar l 



(8) 



where X q = E + m, E is energy of confined quark and m is quark mass. Using the 
substitution p = (Ag^^r for convenience, Eq. (7) reduces to the form 



2 dp 2 



g(p) = eg(p), 



(9) 
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where 



A 2 



i 

i/+2 



(E-m-V ), (10) 
and g(p) is chosen as 

g{p) « pexp (- , (11) 

where x and <i are variation parameters and they are obtained by minimizing e. 
Hence they are solutions of 

(12) 
(13) 



dx ^' 
de 

dd=°- 



Using Eq. (9), Eq. (11) and Eq. (12), e is found to be 

V + 2 N 



e = 



ax 



(14) 



i 

u+2 



n - ( d+l\ o|£l|) 



b = 2^^|^. When Eq. (13) is 



where x = (*) - 

calculated numerically, d is obtained as 1.502 and using this value e is found as 
1.3268 which is very close to the results of the WKB method and 1/iV-expansion. 
Thus, the lowest eigenmodes corresponding to the positive and negative energies 
have the respective explicit forms 



V>A(i Sl/2 )M = 

<t>h{is m ){r) = 



47T 
1 



ig{r)/r 
a ■ rf(r)/r 

a • rf(r)/r 
-ig(r)/r 



Xr, 



Xm,-, 



where the two component spinors Xm and Xm denote xi 
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upper and lower component solutions corresponding to a quark flavor q are 



(15) 
(16) 
and 



. Using Eq. (8) and Eq. (5), the radial parts in the 



9q(l 



f q (r) 



N q ^ UA-[- 



N q d 



r d exp I — 



\r Qq \r 0q J \ \r 0q ^ 

where, N q is normalization constant obtained from the equation 



N *r ( / ' (r)+ ^ 2(r) ) dr=i - 



(17) 
(18) 

(19) 



3 Quark-antiquark momentum distribution 

Knowing the quark-antiquark eigenmodes in the ground-state of meson, it is possible 
to obtain their corresponding momentum distribution amplitude. If G q (p, A, A') is 
the amplitude for finding a bound quark of flavor q in its eigenmode in a 

state of definite momentum p and spin projection A', then it is given by 15 ' 17 

*$\r) = E / d3 p°^ A > ^\jr p u ^ A ') exp f ) • ^ 2 °) 

Where {7 9 (p, A') is the usual free Dirac spinors. Eq. (20) can be easily inverted to 
yield 

G q (p, A, A') = ^Ul(p, A') J <Pr*$\r) exp (-ip- r) . (21) 
Thus, it is found as 

G q (p,\,\') = G q (p)5 xx ,, (22) 

where 



mN q m + E p 



G q (p) = — q —(E p + E q )J —L-J>H{d,z), (23) 

V2dr 0q pi\ q V E p 



where, E p = \Jp 2 + m 2 , E q is the solution of Eq. (10), (3 = {j^j , z = and 

d 2 -3 f°° 2d+i / fv\ d \ 
H(d,z) = -^EE*- J o y~ exp f - (JJ j J 3/2 (yz)dy. (24) 

4 Leptonic decay widths 



Now following Margolis and Mendel 17 one can represent the ground state of a neu- 
tral vector meson such as (p, u, (j)) with a particular spin projection Sy and zero 
momentum as 

\V(0), Sy) = £ j d 3 P G q ( P )C s x ^bl(p, Ai)6j(-p, A 2 ) |0) . (25) 

Here, b q (p, Ai) and b q (— p, A2) operating on the vacuum state are quark and antiquark 
creation operators, respectively. The summation with the flavor coefficient and 
the spin configuration coefficient represents the appropriate SU(6) spin-flavor 

structure of the particular vector meson V with its spin projection Sy and zero 



momentum. The factor is due to the effective color singlet configuration of the 
meson. N(0) represents the overall normalization, which is given by 



1 



DC 



N{0) = j2^fJo d3p \ G «^ 2 - (26) 
Assuming that the main contribution to the leptonic decay process of neutral vector 
mesons such as p, u, (f> comes from single virtual photon creation from the annihi- 
lation of the bound quark-antiquark pair inside of the meson, S'-matrix element in 
configuration space can be written as 15 

S fi = (e~(k l ,5 l )e + (k 2 ,5 2 )\ - ie 2 J d A xid A x 2 (il> e (x2)j fJ 'ipe(x2)D^ u (x2 - x±) 

x J2 e iM x ^ U ^i))\^S v ) (27) 
i 

Where, D fMU (x2 — x\) is the photon propagator, ip e (x) and ip q (x) are the free lep- 
ton and quark fields, respectively. After some standard calculations (details of the 
calculation can be found in Ref. 15), one obtains 

r( V ^e + e-) = —a 2 em M v fl (28) 

Where, jy is known as the leptonic decay constant and in this model it can be 
written as 

2 _ 2{e q ) 2 v I v 
JV ~ TTjTTT?, t ■ t zy ) 



IV 3tt 2 M v J v ' 



Here, 



= (± 



1 1 



Iv = fdpp 2 (2+^G q (p), 

roo 

J v = / dpp 2 \G q (p)\ 2 . (30) 
Jo 

Iv and Jv values can be calculated numerically by computer. 

5 Weak leptonic decay widths 

Here, the weak leptonic decay of charged pseudoscalar mesons such as 7r , 
D ± , Df, B ± and Bf are considered. Assuming that the main contribution to the 



weak leptonic decay processes comes from the single virtual boson creation from the 
annihilation of the bound quark-antiquark pair inside the pseudoscalar meson M, 
the S*-matrix element in configuration space is written as 16 



Sfi = (l(ki,Si)vi(k2,S 2 ) 



-JGf 

7T 



l x^(x) 7 ^(l-7 5 ) Ux) 



E ^^W7m(1-7 5 )^„W|M(0)) 



(31) 



qm,qn 



Where, Gf is the Fermi coupling constant, v qmqn are the CKM matrix elements and 
|M(0)) is given by 16 



|M(0)> = 



N(0) 



^C^ 2 (AiA2) 



x J d 3 p[G 9l (p)G* 2 (-p)] l 6j 1 (p,A 1 )6t 2 (- p ,A 2 )|0), (32) 

where, C^ q2 {\\\2) stands for the appropriate SU(6) spin-flavor coefficients for the 
pseudoscalar meson M. N(0) represents the overall normalization, which is given 
by 



N(0) 



1 



d?p G 51 (p)Gl (-p 



After some standard calculations, which can be found in Ref. 16, one obtains 



T{M^lV l ) = ^\v qiq2 \ 2 M p m 2 (l-|| 
where fu is the weak decay constant, having the form 

ot-2 

JM 2tt 2 M p J m ' 



fit: 



(33) 



(34) 



(35) 



Here M p is the mass of pseudoscalar meson, mi is the mass of lepton and 1m and 
Jm are found as 



f°° r 

/m = J d P p 2 A(p)[G qi (p)G* q2 (-p) 



J M = I dpp 2 G qi (p)G* q2 (-p) 



and 



A(p) 



(E pi + m qi ) {E p2 + m q2 ) - p 2 
[E P1 E P2 (E P1 + m qi ) (E P2 + m 92 )] 2 



(36) 



(37) 



where E pi = ^Jp 2 + m 2 q 



6 Results 



The calculations involve the potential parameters of the model (A,Vo,v) and the 
quark masses (m u = md-,va s ^m c ,m^) . The potential parameters are chosen to be 
A = 0.68GeV, V = -0.3961GeV, v = 0.2. 

The light-quark masses m u = and m s are obtained from to and <fi mesons 
as m u = md = 0.078GeV, m s = 0.3GeV, and from D ± and B ± m c = 1.3GeV, 
m b = 4.81GeV. 

When these parameters are used the masses of D, D s , B, B s , B c ,w and <\> can 
be calculated to be almost the same with their experimental values. Also, with 
the same parameters, radiative decay widths of light and heavy mesons have been 
calculated and the obtained results are close to their experimental values 8 . Using 
this potential parameters and quark masses given above, the calculated results are 
shown in Table I, Table II and Table III. 

7 Conclusion 

In this paper, the independent particle model approach has been used to investigate 
leptonic decay of light vector mesons and weak leptonic decay of light and heavy 
pseudoscalar mesons. The Dirac equation in a power-law potential has been solved 
with variation technique using special type wave function. Leptonic decay widths 
and decay constant fy and weak leptonic decay widths and decay constant Jm have 
been calculated and are compared with the results of other theoretical investigations 
and experiments. The results say to us that the potential used in this study can be 
considered as interaction potential for quarks. Because of the model structure it is 
assumed that quarks do not interact with each other. However, there are spin-spin 
and other hyperfine interactions between quarks. Adding these interactions, good 
results, especially in the calculation of meson mass spectra, could be obtained. Such 
interaction terms can be found in Ref. 24,25. 
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Table I. Leptonic decay widhts and the decay constant fy in keV in comparison 
with the results of other researchers and the experiment. 



Meson 


T{V -> e+e") 


Experiment 23 


Ref. 15 


Ref. 17 


P 


6.37 


6.77±0.32 


6.26(8.1) 


7.8 




0.684 


0.6±0.02 


0.67(0.87) 


0.84 


4> 


1.46 


1.37±0.05 


1.58(1.84) 


1.69 


Meson 


fv 


Experiment 23 


Ref. 15 


Ref. 17 


P 


0.193 


0.2±0.04 


0.19(0.22) 


0.21 


UJ 


0.0624 


0.06±0.01 


0.06(0.07) 


0.07 


4 


0.0813 


0.08±0.01 


0.08(0.09) 


0.07 
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Table II. Decay constants of pseudoscalar mesons in MeV in comparison with 
the results of other model and the experiment. Experimental values are taken from 
Ref. 16. 



Model 


U 


fx 


fD 




fB 


Ib s 




Expt. 16 ' 24 


131.73±0.15 


160.6±1.3 


<219 


137-304 








This work 


131.4 


150.2 


181.1 


235.9 


208.3 


257.3 


377.2 


Ref. 16 


138 


157 


161 


205 


122 


154 


221 


Ref. 14 


100 


153 


149 


160 


96 


111 


141 


Ref. 20 






172 


196 


149 


170 


255 


Ref. 18 


139 


176 


150 


210 


125 


175 


425 


Ref. 19 


178 


182 


148 


166 


98 






Ref. 21 






174±53 


234±72 


105±34 


155±75 




Ref. 22 


141±21 


155±21 


282±28 




183±28 
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Table III. Partial decay widths T (M — > Wf) in MeV and the branching ratio 
B (M — > Ivi) of pseudoscalar mesons in comparison with the experiment. (B (M — ► Ivf) = 
tmP (M — ► /X^)). tm is the mean lifetime of meson M 



Presses 


r(M -► / 


n) 


B (M — > ZFj) 


Expt. 23 B (M — > ZP;) 


TT^ 1 — > U ± V, J 


2.523 10- 


-14 


0.994 


0.999877±0.0000004 


TT^~ — > e^I^e 


3.238 10" 


-18 


1.276 10" 4 


(1.23±0.004) 10~ 4 


K± - ^ 


3.000 10" 


-14 


0.563 


0.6351±0.0018 


^± _^ gi^ 


7.724 10" 


-19 


1.45 10~ 5 


(1.55±0.07) 10~ 5 


£) ± -► ^ 


1.795 10" 


-13 


2.87 10~ 4 


<7.2 10~ 4 


L>± - 


6.240 10- 


-12 


4.41 10~ 3 


4l|;gxl0 -3 


5± ^ ± F At 


1.693 10- 


-16 


4.15 10~ 7 


<2.1 10~ 5 


B± -> ^ ± !7 At 


8.960 10- 


-14 






£)± _> T ±i7 T 


4.720 10" 


-13 


7.540 10~ 4 




£>± -► T ± V T 


6.090 10- 


-11 


4.3 10~ 2 


(7±4) 10~ 2 




3.770 10- 


-14 


9.25 10~ 5 


<0.57 10" 3 


£± - r±I7 T 


2.140 10" 


-11 
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